Minkowski vacua can be metastable 
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Abstract 

We investigate the recent suggestion that a Minkowski vacuum is either absolutely stable, or it has 
a divergent decay rate and thus fails to have a locally Minkowski description. The divergence comes 
from boost integration over momenta of the vacuum bubbles. We point out that a prototypical 
example of false-vacuum decay is pair production in a uniform electric field, so if the argument 
leading to the divergence is correct, it should apply to this case as well. We provide evidence that 
no catastrophic vacuum instability occurs in a constant electric field, indicating that the argument 
cannot be right. Instead, we argue that the boost integration that leads to the divergence is 
unnecessary: when all possible fluctuations of the vacuum bubble are included, the quantum state 
of the bubble is invariant under Lorentz boosts. 



I. INTRODUCTION 



In an interesting recent paper [T], Dvali suggested a simple argument indicating that a 
Minkowski vacuum cannot be metastable. The argument can be summarized as foUows. 

Suppose we have a theory admitting a zero-energy metastable (false) Minkowski vacuum 
which can tunnel to a negative- energy, anti-de Sitter (AdS) true vacuum through bubble nu- 
cleation. The tunneling process conserves energy, so the critical Coleman-De Luccia (CDL) 
bubble has zero energy, with the negative energy in the bubble interior compensated by the 
positive energy of the bubble wall [2j . A bubble whose size is larger than critical and whose 
boundary is at rest will have negative total energy. This state of negative energy may then 
be created from the false vacuum together with a compensating blob of ordinary positive- 
energy particles, perhaps with some separation between the positive and the negative energy 
blobs. 

One might think that the negative-energy blob is like a ghost, that is, a negative-energy 
point-particle. If the energy-momentum vector of the bubble is p^, the positive energy 
blob must have energy-momentum vector — p^, and the Lorentz invariance of the Minkowski 
vacuum implies that the rate for this process per unit spacetime volume can depend only 
on p^. Given a positive plus negative blob at nucleation, we can obtain other possible 
configurations by applying arbitrary Lorentz boosts. This will give configurations with 
arbitrarily high positive and negative energies, adding up to zero. Since does not change 
under boosts, the nucleation rate for all these configurations should be the same. Then the 
total vacuum decay rate, obtained by integrating over the boosts, is infinite. The conclusion 
is that the concept of a metastable Minkowski vacuum is inconsistent; such a vacuum simply 
cannot exist. We shall refer to this as the "boost integration argument." 

As noted in [Ij, a possible caveat is that a metastable Minkowski vacuum can exist only for 
a finite time and thus cannot be Poincare-invariant. This results in a cutoff for the divergent 
boost integral. (This issue will be discussed in more detail in subsequent publications 0.) 
However, according to the standard analysis, the bubble nucleation rate can be arbitrarily 
low, resulting in an arbitrarily high boost cutoff and an arbitrarily large enhancement of the 
vacuum decay rate. This would still call for a drastic revision of the standard approach. In 
what follows we shall use the term "metastable Minkowski vacuum" with the understanding 
that this vacuum is only locally Poincare-invariant. 
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It is important to note that for theories without gravity the argument in [T] does not rely 
in any way on the existence of negative energies. Given a theory with a metastable zero- 
energy vacuum, decaying to a true negative-energy vacuum, we can always add a constant 
to the Lagrangian of the theory, so that both vacua become positive-energy. This of course 
does not change the dynamics, and the argument should still go through, except now the 
energy-momentum vector of the bubble should be understood as the difference between 
the energy momentum vectors of the false vacuum with and without a bubble. If the phase 
space integral in the decay rate was divergent in the original theory, it will still be divergent. 
Hence, if the argument is correct, then all non-gravitational theories with metastable vacua 
should be inconsistent. 

The boost integration argument is reminiscent of the problem encountered in the early 
treatment of vacuum decay [4J by Voloshin, Kobzarev, and Okun (VKO). If we think of 
a CDL bubble as materializing at rest in some frame of reference, this nucleation event 
is not boost-invariant. VKO argued that the nucleation rate should therefore include an 
integration over boosts, resulting in a divergent total rate. (They suggested a cutoff related 
to the radius of the universe.) The key insight by Coleman [5j was that the bubble evolution 
after nucleation is Lorentz invariant. Hence, integration over boosts is unnecessary, since it 
amounts to multiple counting of the same final states. 

Coleman's result is supported by an independent calculation of the vacuum decay rate 
in terms of the imaginary part of the vacuum energy density [6j. In this approach, the issue 
of boost integration does not arise. Another piece of evidence comes from the calculations 
of pair creation rate in a constant electric field [3 ~ ^ process closely analogous to vacuum 
decay. Starting with the work of Heisenberg and Euler [8j and Schwinger [9j , pair creation 
has been studied by a variety of methods. Results obtained using different methods are in 
full agreement and do not include any divergent boost integrals. 

Here we are going to argue that Dvali's conundrum can be resolved in a similar way. 
An expanding bubble can emit particles, so the asymptotic states at future infinity do 
include positive-energy particles with a compensating negative-energy bubble. In fact, the 
total number of emitted particles and their total energy are going to be infinite in this 
limit. Suppose the bubble were formed together with a photon, as in the boost integration 
argument. This process cannot be distinguished from bubble formation with a subsequent 
emission of the photon at about the same time. Now, if we apply a Lorentz boost, the 
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emission point of the photon wiU move to the future (or to the past). Large Lorentz boosts 
wiU move it to the very distant future. Integration over boosts would thus account for aU 
possible photon emission points. Now, the probability that the bubble will emit a single 
photon in its entire history is zero. As we said, the total number of emitted photons is infinite, 
and all final states having nonzero probability correspond to a certain average number of 
photons emitted per unit area of the bubble wall per unit proper time. Such states remain 
unchanged under Lorentz boosts. 

Another way of looking at this is to note that the bubble wall has the geometry of a 
(2+l)-dimensional de Sitter space. Any quantum field interacting with the field(s) of the 
bubble will be in a de Sitter invariant quantum state around the wall [lOj. This quantum 
state accounts for all particles that will ever be emitted by the wall. Since the de Sitter 
space is characterized by an intrinsic temperature, things more substantial than particles 
will sometimes also pop out. There will be blobs of matter, occasional Boltzmann brains, 
etc. This will be happening all along the domain wall in a de Sitter invariant fashion. 
Now, Lorentz boosts in the bulk spacetime correspond to de Sitter transformations on the 
worldvolume, and since the quantum state is de Sitter invariant, it does not change under 
boosts. Hence, including configurations of negative and positive energies does not break the 
Lorentz invariance of the final state, but is necessarily part of the quantum description of a 
Lorentz invariant bubble. 

To substantiate this view, we shall investigate how the boost integration argument plays 
out in the case of pair production in electric field. We shall first focus on the (1+1)- 
dimensional case, which is a prototypical example of vacuum decay. Starting with a brief 
overview in the next section, we shall analyze some simple massless (1+1)D models in Section 



III, Massive particle production in (3+l)D is discussed in Section MVl and our conclusions 



are summarized in Section [Vl 



II. PAIR PRODUCTION IN AN ELECTRIC FIELD 

A constant electric field in (1+1) dimensions, = Ecj^jy with E = const, is invariant 
under Poincare transformations. Its energy-momentum tensor has the vacuum form, = 
If the field is coupled to particles of electric charge e and mass m, particle-antiparticle 
pairs will be spontaneously produced with a critical separation dc = 2m /eE^ so that the rest 
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energy of the particles, 2m is compensated by the potential energy, —eEdc- (For definiteness 
we assume that E and e are both positive.) The electric field in the space between the pair 
is reduced to the value (£" — e), and the vacuum energy is reduced accordingly. The created 
particles play the role of bubble walls in this model. The rate of pair creation per unit length 
is (see, e.g., Ref. [llj and references therein)-^ 

eE f 7rm^\ 

If the particles are coupled to themselves or to another particle species, then pairs sep- 
arated by a distance d > dc can be produced. Such pairs have negative energy, and their 
formation must be accompanied by production of compensating positive-energy particles. 
It is clear that the boost integration argument is precisely concerned with this situation. If 
the argument is correct, the pair production rate should diverge when particle interactions 
are included. We shall now examine whether or not such a divergence actually occurs. 



III. MASSLESS PAIRS 



Massless QED in (1+1) dimensions, also known as the Schwinger model, is described by 
the action 



S= d^x 



It is equivalent to the bosonic theory [T3] 

"1 



S= d'x 



where g e/y^. 

In terms of the boson, the electric current is given by 



(2) 



(3) 



(4) 



The number density of charge carriers n is related to the electric current j = —gcj) by 
j = 2e n. It follows that the pair production rate is given by 



2edt 20?' 



(5) 



^ This expression is exact in external field (corresponding to the hmit e <^ E) . Backreaction can easily be 
taken into account in the instanton approximation (see, e.g., [12 ). In this case, the result is given by (1), 
but with E replaced hy E = E — e/2. The instanton approximation is valid in the hmit ^ eE. 
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The field equation for (j) is 

□0 = -gE, (6) 
where E is the electric field. For a constant external field, (j) develops an expectation value 

(0) = ut) = -IgEf, (7) 



and Eq. g sives 



e 



As one might expect, there is no exponential suppression of pair production for massless 
fermions. 

We note that the massless pair production rate ([s]) is in agreement with the zero-mass 
limit of Eq. Q. (For a detailed discussion of this limit, see [llj.) It was first derived by 
Witten [14j in the context of superconducting cosmic strings. 

In this paper we will not be interested in the effect of back-reaction of the created pairs 
on the electric field. However, if needed, this effect can be easily taken into account pTl [T5| . 
From Maxwell's equation 

E = -j = g4>, (9) 



and using (j6]) we have 



and 



E = Eo + g4> (10) 



D(t> + g^(t>=-gEo, (11) 



where £"0 = const. The solution of this equation for a spatially homogeneous (/) with the 
initial condition = = at t = 0, is 

Mt) = -9~'Eo[l-cos{gt)l (12) 



and Eq. (10) gives 

E = Eoco8{gt). (13) 

This shows that the screening of the electric field by the produced pairs occurs on a timescale 
t ^ g~^. At smaller times t <C g~^^ Eq. (12) is well approximated by Eq. ([t]). 



The above results are valid to all orders in e, and show that catastrophic vacuum decay 
does not occur in the Schwinger model. One may be concerned that this model is somewhat 
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special. Note that there are no propagating photons in (1+1) dimensions. Hence, if we 
neglect the dynamics of the electric field, the model does not include any particles that 
could compensate for the increased negative energy of the pairs. On the other hand, if we 
do include a dynamical electric field, the model is still equivalent to a free bosonic theory, 
and hence it may be considered non- generic. In order to account for more generic situations 
with propagating particles in final states, we now consider more complicated interacting 
theories. 



A. Interactions 

The Schwinger model ^ can be modified by adding self-interaction terms for the field 
ip. We first consider a simple interaction of the form 

>C/ = A j,f, (14) 

where A is a coupling constant. Apart from the Maxwell field, the action ^ with this 
interaction term corresponds to the massless Thirring model. Naively, the four fermion 
interaction could lead to pair creation processes which are not boost invariant, so one might 
expect catastrophic vacuum decay once an electric field is applied. 

This, however, does not happen. In the bosonized form of the model, the interaction 



term (14) becomes 

Ci = X9'{d<pf, (15) 

and amounts to a finite renormalization of the kinetic term for (j). Thus the model is 
equivalent to the bosonic Schwinger model, with the replacement g ^ = {1 + 2\g^)~^/'^g. 
If we apply a constant electric field, the electric current will grow at a constant rate, as 
before: 

dj ^ e eE 

dt 27r(l + 2A^2)- V ; 

Clearly, no catastrophic vacuum decay occurs in this model. 

Turning now to more generic interactions, we introduce an additional scalar field x, which 
is coupled to ip via 

Cj = H{x)e'^d,j,, (17) 

where the function H{x) is assumed to be 0{x^) at small x a^id is otherwise arbitrary. This 
form of interaction is chosen so that the model can be rewritten in an equivalent bosonic 



form. We could also choose a coupling G{x)jiij^^ but since the current has an expectation 
value linearly growing with time (at least to the lowest order in the interaction), this would 
result in a time-dependent, growing mass for x- 

Using Eq. Q to express the current in terms of the bosonic field 0, we have 

e^'^d.j, = 2^00, (18) 



so the bosonic action takes the form 

1 



S= d'x 



my + [dxY] - V{x) + H{x)n<l) - gE(/> 



(19) 



where we have also added a self-interaction potential for x- 

To lowest order in the interactions, the expectation values of (/) and x are {(/)) = 00 (t) and 
ix) = 0, with (j)o{t) from Eq. ([t]). Introducing = — ^q, the action transforms into 

"1, 



S 



I 



d X 



+ {dxY] - V{x) + H{x){U4> - gE) 



(20) 



For a constant electric field, the term proportional to gE amounts to a finite renormalization 
of the mass and self-couplings of x- 

We see that even though the classical solution 00 (t) is growing with time, the model 



(20) describing the quantum theory on this background has a time-independent Lagrangian 
and shows no signs of instability. Once again, it seems clear that nothing catastrophic can 
happen to the vacuum in this model. 



IV. MASSIVE QED 

We finally consider pair production in massive QED, with photons playing the role of 
the compensating positive-energy particles. Since there are no photons in 2D, we consider 
the 4-dimensional case. A constant, spatially homogeneous electric field is invariant under 
longitudinal boosts, so according to the boost integration argument, the rate for processes 
like 

Vacuum e^e~7. Vacuum e^e~77, (21) 

etc., should include integration over such boosts and should therefore be divergent. 

The total vacuum decay rate, including all particle production processes, is related to the 
imaginary part of the vacuum energy density, 

Tvac = 2Impvac, (22) 
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(a) 



(b) 



(c) 




(d) 



FIG. 1: Diagrams contributing to the decay of an electric field. The solid lines represent the 
electron propagator in the presence of the electric field, and the wavy lines represent the photon 
propagator. 

and can be found by evaluating the diagrams shown in Fig. 1. Solid lines in these diagrams 
stand for electron propagators in a constant external field, and wavy lines represent photon 
exchange. The one-loop diagram in Fig. 1(a), which does not include any virtual photon 
lines, gives the Heisenberg - Euler - Schwinger (HES) result. 



r(i) 



{eEf 

47^3 



exp 



Tim 



(23) 



The higher-loop diagrams of Figs. l(b)-(c) account for the processes (21), where pair pro- 



duction is accompanied by photons, as well as for the radiative corrections to the basic 
pair production process. Vacuum e+e~. The diagram in Fig. 1(d) represents correlated 
formation of two pairs and radiative corrections to the photon propagator. If the boost 
integration argument is correct, the multi-loop diagrams should diverge. 

It should be noted that the pair production rate F is generally different from the vacuum 



decay rate Fvac defined by Eq. (22). At the one-loop level, the exact result for F is given by 



Eq. (23), while F^ac is |9] 



F = 



[eEf 



exp 



rnim 
eE 



(24) 



For an illuminating discussion of this distinction, see, e.g., Ref. [TT].^ Here, we will be 
interested in the weak field limit, 

E < mVe, (25) 

when r ^ Fvac to very good accuracy. Hence, we will not distinguish between Fyac and F in 
what follows. 

Higher-loop corrections to the HES formula ([23]) have been extensively studied in the 
literature. Affleck, Alvarez and Manton (AAM) [16] derived a remarkably simple formula, 
which includes contributions of all multi-loop diagrams containing a single electron loop, 

r.--^<g!exp(-^4). (26) 

It was argued by AAM that diagrams with more than one electron loop are subdominant in 



the weak field limit (25); then Eq. (26) is the full multi-loop pair creation rate in this limit. 
Comparing it to (23), we see that the correction term is definitely finite and small, as one 
might naively expect. 

The AAM derivation was not rigorous, as it relied on the steepest descent evaluation 
of the path integral for the QED effective action and on a somewhat cavalier treatment of 
mass renormalization. However, direct calculations p!7l - [T9] of the two-loop contribution in 



Fig. 1(b) also yielded a finite result, which is in agreement with (26). 

It may be instructive to consider how the infinite boost integration plays out in the QED 
context. In the semiclassical approximation, electrons (and positrons) move along hyperbolic 
trajectories with a constant proper acceleration, a = eE/m. The characteristic energy of the 



emitted photons in the rest frame of the electron is ~ a, and in the weak field limit (25) 
we have <C m. This indicates that photon emission has negligible effect on the motion of 
the electrons, and one can simply consider radiation from a charged particle moving along 
a fixed classical trajectory. 

For a charge in hyperbolic motion, photons will be emitted with the same characteristic 
energy in the rest frame of the charge, but will be boosted to arbitrarily high energies 
in the observer's frame. However, the transverse momentum of the photons, in the plane 
perpendicular to the electric field, will not be boosted. Hence, the appropriate quantity to 



^ The distinction is particularly important in the massless (l+l)-dimensional case, when Fvac becomes 
formally divergent, while the pair production rate (p| remains finite [llj . 
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calculate is the probability of emitting a photon with a given transverse momentum p^. 
This calculation has been done, in the lowest order of perturbation theory, in Ref. [20] for 
photons and in Ref. [21 J for the emission of massless spin-zero particles by a scalar source. 
The resulting probability is proportional to the divergent integral over the particle's proper 
time (which can also be thought of as an integral over boosts, since a boost causes a shift 
in the proper time of emission). How do we interpret this divergence? 

Of course, the probability of photon emission cannot be greater than one, so the diver- 
gence indicates a breakdown of the perturbation theory. The origin of the problem can be 
understood if we consider a situation where the electric field is turned on only for a finite 
period of time T. We shall assume this period to be short enough, so that the probability 
of emitting more than one photon in time T is small. (Note that we can treat the electric 
charge e as a free parameter. For small values of e, photons will be emitted very rarely, so 
the time T can be very long.) In this case there will be no divergence. The probability of 
photon emission will be proportional to the proper time r spent by the charge in the electric 
field. This simply refiects the fact that there is a fixed photon emission probability per unit 
proper time. 

What happens when we increase the time T? For very large values of T, the probability 
of emitting a single photon becomes very small. Photon emission events along the particle's 
trajectory can be regarded as independent, so the probability for emitting n photons will be 
given by the Poisson distribution 

iy=.n 

Pn = -^e-^ (27) 

ni 

where n is the average number of the emitted photons (which is proportional to r). In the 
limit T ^ oc, we have n ^ oc, so the probability of emitting any finite number of photons 
goes to zero. On the other hand, n is proportional to e^, so a formal perturbative expansion 
gives a divergent result for the emission of a single quantum, p^^^^^ = n + 0{e^).'^ 

^ The situation here is similar to the divergence encountered in emission of soft photons, known as the 
infrared catastrophe in QED; see, e.g., [22] . 
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V. CONCLUSIONS 



In this paper we have analyzed Dvah's boost integration argument indicating that a 
Minkowski vacuum can either be absolutely stable or cannot exist at all. Any decay channel 
allowing negative-energy bubbles would supposedly result in a divergent vacuum decay rate. 
This would imply good news for inhabitants of a Minkowski vacuum: if your vacuum survived 
for any finite time, it is guaranteed to endure forever. Since our vacuum is close to Minkowski, 
by continuity we should then expect it to be extremely long lived. 

Here, we argued that if the boost integration argument is correct, it should also apply to 
pair production of charged particles in a constant electric field. The pair production rate 



should then become infinite when processes like (21), where the pairs are produced together 
with some other particles, are included. 

We do not know of any calculations that would exhibit such a divergence. On the other 
hand, here we have presented evidence that such divergences do not occur in a class of 
massless (l+l)-dimensional theories and in (3+l)-dimensional (massive) QED. 

The divergence found in the boost integration argument arises due to the integration over 
the momenta of the vacuum decay products. A perturbative calculation of the probability 
for nucleating a negative energy bubble accompanied by a single positive energy blob of 
matter in the final state contains a divergent integral over the momentum of the blob. 

However, as we discussed in Section IV, this divergence only indicates a breakdown of per- 
turbation theory. Physically, the reason is the following. Fluctuations such as the emission 
of a blob by an expanding bubble occur with a constant probability per unit wall area per 
unit proper time. The boost integral in the naive perturbative calculation only changes the 
value of proper time on the worldsheet at which a blob is emitted. If the bubble expands for 
a finite time, i.e., short enough that the emission of one blob is unlikely, then perturbation 
theory applies, and the boost integral will correctly account for the fact that the probability 
of emission is linear in the proper time interval available. However, it is clear that this inte- 
gral is unrelated to the probability of nucleating the bubble itself. As the bubble expands, 
the probability of having the bubble with the blob increases, and the probability of having 
the bubble without the blob decreases. Once the probability for having more than one blob 
is significant, perturbation theory has broken down. 

In a Lorentz invariant situation, where the bubble expands for an infinite amount of 
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time, the correct answer for the probabihty of having just one blob (or any finite number of 
them) in the final state is actually zero. In other words, the final state consists of a Lorentz 
invariant distribution of infinitely many blobs, generated as the bubble expands into the 
false vacuum. But of course this does not imply an infinite rate of bubble nucleation. 

A metastable Minkowski vacuum cannot exist forever and thus cannot be truly Poincare 
invariant. However, such vacua can be formed as bubbles in the infiationary multiverse and 
will decay by producing AdS bubbles. Our analysis indicates that the corresponding decay 
rate can be calculated using the standard Coleman-de Luccia method and does not involve 
any divergent boost integrations.^ 
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